Long-Range Interaction of Spin-Qubits via Ferromagnets 
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We propose a mechanism of long-range coherent coupling between spins coupled to a ferromagnet 
by exchnage or dipolar coupling. An effective two-spin interaction Ifamiltonian is derived and the 
coupling strength is estimated. We also discuss mechanisms of decoherence and consider possibil- 
ities for gate control of the interaction between neighboring spin-qubits. The resulting quantum 
computing architecture retains all the single qubit gates and measurement aspects of earlier ap- 
proaches, but allows qubit spacing at distances of order 1 /im for two-qubit gates, achievable with 
current semiconductor technology. The clock speed depends strongly on the dimensionality of the 
ferromagnet and is between MHz and GHz. 

PACS numbers: 73.20.Dx, 71.70.Ej, 03.67.Lx, 76.60.-k 



Introduction. Quantum coherence and entanglement 
lie at the heart of quantum information processing. One 
of the basic requirements for implementing quantum 
computing is to generate, control, and, measure entan- 
glement in a given quantum system. This is a rather 
challenging task, as it requires to overcome several ob- 
stacles, the most important one being decoherence pro- 
cesses. These negative effects have their origin in the 
unavoidable coupling of the quantum systems to the en- 
vironment they are residing in. 

An unwritten rule in the search for a good system to 
encode qubits is the smaller the system the more coher- 
ence, or, more precisely, the fewer degrees of freedom the 
weaker the coupling to the environment. Simultaneously, 
one needs to be able to coherently manipulate the indi- 
vidual quantum objects, which is more efficient for larger 
systems. This immediately forces us to compromise be- 
tween manipulation and decoherence requirements. 

Following this rule, one of the most successful candi- 
dates for encoding a qubit is an electron spin localized 
in a semiconductor quantum dot, gate-defined or self- 
assembled [H. Indeed, this natural two- level system ful- 
fills most of the prerequisites for a good qubit; it is very 
long-lived (relaxation time Ti ~ Is Q and decoherence 
time T2 '-^ 300/Lts it can be controlled efficiently by 
both electric and magnetic fields and, eventually, 

may be scaled into a large network. However, one of the 
central challenges for this scaling remains the qubit-qubit 
coupling, in particular, via distances that give enough 
space to accommodate gates and control lines that are 
needed for engineering large networks. 

Addressing this issue, we propose in this work a setup 
to couple two spin qubits separated by a relatively large 
distance on the order of micrometers, see Fig.[T^. In con- 
trast to other schemes [7-ll|, the coupling is mediated 
via a ferromagnet to which the spin qubits are coupled by 
exchange or dipolar interaction. The local exchange cou- 
pling between the ferromagnet and the spin qubits can 
be realized by etching a trench close to the spin qubit 
and placing a ferromagnet therein, see Fig. [TJd. Alterna- 
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FIG. 1. The schematics of the two-dimensional quantum com- 
puter architecture. The left panel (a), shows the top view of 
the architecture proposed. The orange parts represent fer- 
romagnet, while the blue circles represent localized electron 
whose spin (green) is used to define a qubit. The right panel 
(b) depicts the side view of the setup. The ferromagnet can 
be placed in vicinity of the qubit by etching a trench and thus 
also removing part of 2DEG (black line). 



tively, the ferromagnet can be placed on top of the qubit 
structure and coupled to the qubit spin via dipolc-dipole 
interaction. The switching on and off the qubit-qubit 
interaction is achieved either by controlling the tunnel 
coupling to the ferromagnet by electrical means or by 
tuning qubits off resonance (see below). The resulting 
system is thus realizable with state-of-the-art semicon- 
ductor technologies. 

Let us first give an intuitive picture of the qubit-qubit 
coupling, before we proceed with the quantitative anal- 
ysis. The coupling between two distant qubits is medi- 
ated via a coupler system. The relevant quantity of this 
coupler is its spin-spin susceptibility — in order to have a 
long-range coupling, a slowly spatially decaying suscep- 
tibility is required. Thus, the dimensionality of the cou- 
pler plays an important role since, in general, it strongly 
influences the spatial decay of the susceptibility, which 
can be anticipated from purely geometric considerations. 
Furthermore, since we are interested in coherent inter- 
actions between the qubits, care must be taken that the 
coupler does not introduce detrimental decoherence — the 
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quantity containing information about this is the noise 
power spectrum S{uj) of the coupler, which, in turn, is 
related to the imaginary part of the coupler susceptibility 
via the fluctuation-dissipation theorem. 

When the qubits are separated by a distance of l^m, 
the coupling strength of qubit-qubit interaction reaches 
values of 10~^°eV for a 3D tunnel coupled ferromagnet. 
Furthermore, we show that the tunnel and dipolar cou- 
plings to ferromagnet do not cause additional decoher- 
ence problems since the decoherence times due to it can 
be tailored to be as long as several minutes. Much better 
qubit-qubit coupling strengths, on the order of lO^^cV, 
can be obtained if a quasi- ID ferromagnet (magnetic 
semiconductor) is used as a coupler. Finally, if the cou- 
pling to the ferromagnet is of dipolar origin, a coupling 
strength of 10~*eV can be reached. 

Exchange coupling. We denote by Sr the spins (of 
size S) of the ferromagnet at site r on a cubic lattice and 
(Ti stands for the spin-1/2 qubit spins. The Hamiltonian 
we consider is of the following form 



by employing a perturbative Schrieffer- Wolff transforma- 



H = Hp 



(1) 



where A is the exchange coupling constant between the 
qubit spins and the ferromagnet (for the dipolar cou- 
pling, sec below). The Hamiltonian of the ferromagnet. 
Hp — ~Jj2{rr')'^r ■ With J > 0, is thc thrcc- 
dimensional (3D) Heisenberg model with the sum re- 
stricted to nearest-neighbor sites (r, r'). Thc ferromagnet 
is assumed to be below thc Curie temperature with thc 
magnetization pointing along the z-direction. The qubit 
Hamiltonian is assumed to be without splitting initially, 
that is Ha ~ 0. Nevertheless, since the ferromagnet is 
in the ordered phase, there exists a first order effect due 
to coupling to the ferromagnet H^ = AJ2i'^i i^rt) ■ 
will become clear later, such a splitting is undesirable 
if one is interested in coherent interaction — we remedy 
this by coupling thc spins to another ferromagnet, al- 
beit with anti-parallel magnetization. Since we allow for 
some misalignment between orientation of the magneti- 
zation of the two fcrromagncts, thc final Hamiltonian for 
thc qubits in thc spin space after correctly taking into 
account thc first order corrections due to coupling to thc 
ferromagnet reads 



Hr. 



(2) 



The splitting in the x-direction of the qubit (or equiva- 
lently along the y-direction) is beneficial since it reduces 
decoherence due to longitudinal noise of the ferromagnet: 
The effect of such noise spectrum can significantly infiu- 
ence decoherence times for the case of no splitting of thc 
qubit because thc longitudinal noise is gapless. 

Coherent coupling. We proceed with thc derivation of 
an effective two-spin interaction Hamiltonian for A <^ J 



tion 



12l | up to the second order (see Appendix) 



H, 



H„ 



-X±(A)(2a^4-aJaf-a^<72-), (3) 



where we introduced the notation x±i^) = X±{^tL) 
[L = |r2 — ri|) and x±{u), r) is the transverse real space 
spin susceptibility of thc ferromagnet. Note that we have 
neglected Xj^(~A) and x^±^{^) in comparison to Xj^(A), 
as well as the longitudinal susceptibility X|| since it is 
smaller by factor of 1/ S compared to the transverse one 
(sec Appendix). Thc real space transverse susceptibility 
of thc 3D ferromagnet is given by (see Appendix) 



5- e-rV^ 

15 ~ 



00 < Ap, 



(4) 



where Ai? is the gap induced via applied external mag- 
netic field or due to internal anisotropy of the ferromag- 
net, Ip ~ \j Ap-ijj ^'^'^ ^ ^ 2JS. In what follows, 
we assume that the external gap is always larger than 
the qubit splitting, A < A^, as this ensures that the 
transverse noise is not contributing to decoherence since 
transverse noise is related to the vanishing imaginary part 
of the transverse susceptibility, X-l(w)" = (w < Ap). 
The spatial dependence of the effective two spin coupling 
given by Eq. Q is of Yukawa type due to presence of 
the external gap. If we assume a realistic tunnel cou- 
pling to the ferromagnet of lOO^eV Q, a Curie temper- 
ature of lOOK and a gap of Ai? = 100/LteV, and the qubit 
splitting close to the resonance A/? — A = 3 x 10"^ fieV 
(corresponding to a magnetic field of about B = 60^T) 
we obtain for the qubit-qubit coupling strength a value 
on the order of 10~"'^°eV for a lattice constant of about 
4A. This coupling strength gives rise to the operation 
times on the order of 5/is — significantly below the relax- 
ation and decoherence times of the spin qubit, Ti = Is 
and = 300ns [| respectively. Furthermore, the error 
threshold — defined as the ratio between the two-qubit 
gate operation time to the decoherence time — we obtain 
with such an operation time is about 10~^, which is good 
enough for implementing the surface code error correc- 
tion [ill. 

The dimensionality of the ferromagnet plays an im- 
portant role — if we assume lOnm width of the trench 
where the ferromagnet is placed, then, for energies be- 
low O.lmeV, the ferromagnet behaves as quasi one- 
dimensional (ID). In this case we obtain (see Appendix) 



r/lf 



LU < Al 



(5) 



wherefrom it is evident that at distances r < Ip the 
susceptibility of a quasi-lD ferromagnet is practically 
constant in contrast to the 3D case, where a 1/r de- 
cay is obtained, see Eq. Additionally, we require 
Ip < D/{AS) — 2J/A for thc perturbation theory to 
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be valid. Thus, for the same parameters as above, but 
without the need to tune very close to the resonance (wc 
set herein A^- — A = 0.5/ieV, corresponding to about 
B = 5mT) a coupling strength of 4 x lO'^^eV is ob- 
tained. 

For ID case there is yet another rather promising 
possibility — to use magnetic semiconductors These 
materials are characterized by a particularly low Curie 
temperature of 30K or below Q, and the distance be- 
tween the ions that are magnetically ordered via RKKY 
interaction is about 10 — lOOnm. Such a large lattice con- 
stant is very beneficial for the long range coupling — if we 
take the lattice constant to be lOnm, the coupling to the 
ferromagnet A = 15/zeV and the qubit splitting close to 
resonance (Ap — A = 0.5/ieV, corresponding to about 
B = 5?T!T), the qubit-qubit coupling becomes of the or- 
der of lyueV. Such a coupling strength in turn leads to 
an error threshold on the order of 10^^. Therefore, even 
the standard error correction protocol can be used in this 
case. 

Finally we note that the bottleneck for the opera- 
tion time is given by the magnon velocity. The typical 
magnon energy is given by A, and thus the magnon ve- 
locity can be estimated as vm = V JSAa/h, leading to 
times of 10ns (Ins) for magnons to go over distances of 
Ifim in ferromagnets (magnetic semiconductors). 

Dipolar coupling. The exchange coupling to the fer- 
romagnet that we considered in the previous section is 
realistic, albeit experimentally challenging. This moti- 
vates us to consider a more feasible setup wherein the 
ferromagnetic coupler is simply placed on top of the spin- 
qubit structure Q . Since there is no overlap between the 
qubit wave function and the ferromagnet spins in this 
setup, the qubit-fcrromagnet coupling is given by dipole- 
dipolc interaction. Thus, the total Hamiltonian reads 
H = Hp + Her + Hi, where the interaction between the 
qubit and the ferromagnet spin. Hi, reads 

^/ = L^^(^--^r j, (6) 

where /x is the magnetic moment of the ferromagnet spins. 
In order to obtain a sizable coupling and also to sim- 
plify the analysis we assume a 'dogbone' shape of the 
ferromagnetic coupler as in Ref. i—two ferromagnetic 
discs connected via a quasi- ID ferromagnetic line. After 
performing a SW transformation wc obtain the effective 
Hamiltonian 

i/eff = i/. + ^xl°(A, mo\al - alal - alal) , (7) 

where again we retained only the on-resonance suscepti- 
bility, i.e., we neglected xi°(— A) and xi^(O) in compar- 
ison to xl°(A). Here, B is given by (see Appendix) 



where ^2 — di is the thickness of the discs, di is the dis- 
tance in z-direction from the qubit spin to the bottom of 
the adjacent disc, and Rq is the disc radius. If we use the 
same parameters as before for the quasi-lD ferromagnet 
with lattice constant of 4A and Ap — A = 10~^/xeV — a 
coupling strength of 10~^eV is obtained. Here, we as- 
sumed ^1 = lO^iB, di ~ 50mn, d2 — di = 20nm, and 
i?o = 50nm. 

The qubit-qubit interaction in case of the dipolar cou- 
pling to the ferromagnet can be efficiently switched on 
(off) by tunning the qubits on- (off-) resonance. Indeed, 
the susceptibility depends sensitively on the difference 
Aj;- — A, sec Eq. ([5]), it is enough to detune by less than 
1% to practically switch off the coupling. This detun- 
ing can be for example achieved by changing the qubit 
splitting by a spatially dependent g-factor which can be 
controlled by electrical means. Such a switching scheme 
will not decouple the ferromagnet as a source of decoher- 
ence. Nevertheless, as we will demonstrate below, when 
the transverse noise is gapped the decoherence caused by 
the ferromagnet becomes negligible, and thus there is no 
need to decouple the qubit from the ferromagnet during 
the ojf state. 

Two-qubit gates. Since the Hamiltonian of Eq. ([3]) is 
entangling, it can be used to implement two-qubit gates. 
Here wc consider the iSWAP, CPE, and CNOT gates, 
the last one being widely used in schemes for quantum 
computation [l^. The Hamiltonian for two single-QD 
qubits interacting via the ferromagnet is the sum of the 
qubit-qubit interaction and the Zeeman terms, given in 
Eq. ([3]) . The strength of the latter in comparison to the 
former allows us to approximate the Hamiltonian (see 
Appendix) 

H^ff = Ji2(a^f + arV2^) + i/^, (9) 

with the notation J12 = ^^X-l(A)/4 and H^ given by 
Eq. ([2]). After this approximation, the qubit-qubit inter- 
action and Zeeman terms commute. The iSWAP gate, 
C^iSWAPj may then be realized with the following sequence 

c/iswAP = nin2e'"''e""''«'nin2 , (lo) 

where t = 37r/(4Ji2) and H denotes the single-qubit 
Hadamard rotation. More details on the construction 
can be found in the Appendix. Previously, the sequence 
for the universal CPF gate was constructed [l^ for the 
same Hamiltonian as the projected one given in Eq. ([9|), 
albeit without the Zeeman terms which commute with 
the Hamiltonian. With the iSWAP the CNOT gate can 
be constructed in a standard way 

C/cNOT = e^'T'^Je^f '^2e*f '^2[/iswApe"'^'"i t/iswApe'^'^= . 

(11) 

Since H^f^ is only an approximation of the total Hamil- 
tonian, these sequences will yield approximate iSWAP 
and CNOT. Their success can be characterized by the 
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fidelity, as defined in the Appendix. For realistic pa- 
rameters, with the Zeeman terms two order of magni- 
tude stronger than the qubit-qubit coupling, the above 
sequence yields fidelities for the CNOT gate of 99.976%. 
This is well above the fidelity of 99.17%, corresponding to 
the threshold for noisy CNOTs in the surface code [18| . 

Decoherence. Next, we study the dynamics of a single 
spin qubit coupled to the ferromagnet, but neglecting cor- 
related noise (see Appendix). The qubit dynamics can be 
described by the reduced density matrix pR{t) = tr (yo(t)), 
where p is the density matrix for the total system and the 
trace is taken over the ferromagnet. Here, pn{t) satisfies 
the equation [l^ 



PR = -i[Ha,Pvi] + 1^PB 



(12) 



with 2? given by the following expression in the Born- 
Markov approximation 



T^pRit) 



(A)(pR(t)-a^PR(t)a; 



+ 



(13) 



+ A^q{A) {[a;,pR]+ - laypRa; + la^pRal) 

where S\\{uj) is longitudinal noise power of the fer- 
romagnet, defined as Fourier transform of S\\{t) = 
([S'^(t),5^(0)]+), while C|'|'(a;) is defined as imag- 
inary part of the correlation function Cy (a;) = 
dtsin{u;t){S'it)S'{0)) = dtsm(ujt)Cii{t). In writ- 



ing Eq. (|T3|) we assumed that the external gap of the fer- 
romagnet, Ai?, is bigger than the splitting of the qubit 
(i.e. A < Ap). Therefore, the transverse noise of the 
ferromagnet is not entering in the Born approximation 
since it is gapped. The 3D ferromagnet behaves as sub- 
ohmic bath Jy (a;) = ujs with s = 1/2, uis = T and 

Us = aXyD^e^^^^, where a = ^-''f'f<^^ = 0.21 (see 
Appendix). From Eg. we obtain the relaxation and 
decoherence times [l9| for the 3D ferromagnet under the 
assumption Ap ^ T 



t; 



A^5|r(A) 



'^'^^coth(/3A/2). (14) 



The decoherence time becomes T2 = 2Ti. The above 
formula allows us to estimate the additional decoherence 
and relaxation due to the ferromagnet. If we assume tem- 
peratures below O.IK, we get decoherence times due to 
ferromagnet of about an hour. Such a long decoherence 
time is obtained mainly because we choose parameters in 
such a way that the transverse susceptibility is gapped. 
The remaining longitudinal noise of the ferromagnet van- 
ishes at zero temperature and by choosing reasonably low 
temperatures (compared to Ap) this noise source can be 
significantly reduced. Furthermore, for quasi-lD ferro- 
magnets, we obtain divergent (as a; — > 0) spectral den- 
sity with s = -1/2, ujs = T and w,, = ^e"'^'^^. The 
relaxation rate is now given by 
7^2 



A^S}P{A) 



e^''^^ coth(/3A/2) , (15) 



where 7 is a numerical factor of order 1. As before, the 
decoherence time reads T2 = 2 Ti . With the same tem- 
perature and the external gap as for the 3D case, the 
decoherence time is about milliseconds. 

Similarly, the relaxation time for dipolar coupling 
reads 



2a6 

X 5'j|(A,r"-r') 



dr'dr" I ^ iz'^ 



«2 ' 



disc ' ' 



(16) 



Instead of explicitly calculating the above integral we can 
readily give a lower bound for the relaxation (and deco- 
herence) time using a rather general inequality, 5j| . Thus 
the lower bound for the relaxation time is given 



Tf 1 < B^SF{A), 



(17) 



where B is given in Eq. ^ and 5'||°(A) in Eq. dH]). This 
leads to the lower bound of the relaxation time which is 
on the order of minutes. 

We note that the effect of the noise is, by and large, 
more pronounced for A = and this was the main mo- 
tivation for us to introduce A 7^ 0. In order to clarify 
the behavior of the system when A = 0, we study in the 
Appendix an exactly solvable model which demonstrates 
the limitation of the Born approximation. 

Conclusions. The clock speed of the quantum com- 
puter architecture just described appears to be a fraction 
of MHz (3D ferromagnet) and IGHz (quasi ID magnetic 
semiconductor) and should be compared with the time 
scales for relaxation and decoherence. The leading mech- 
anism for these at low temperatures is through interac- 
tion with nuclear spins for, say, spin qubits in GaAs dots. 
Typically, the decoherence time Tj* of order 300^s 0] . 
In our case, there is no significant additional decoherence 
mechanisms caused by tunnel or dipolar coupling to the 
ferromagnet. 

The standard quantum error correction protocols re- 
quire the quality factor, equal the ratio of the gate- 
function clock time to decoherence time, not to exceed 
10~^. Our estimates indicate that this is not the case 
for the present system when a 3D ferromagnet is used. 
However, there exist error correction schemes with sig- 
nificantly higher error threshold, such as surface code 
with error threshold of 10^^ [l3|- Since the architecture 
proposed herein satisfies all the criteria for implementing 
the surface code error correction, we conclude that our 
proposal satisfies the criteria for scalability of a quan- 
tum computer. Additionally, when quasi-lD magnetic 
semiconductors are used as couplers, even the standard 
quantum error correction protocols can be used. 
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HOLSTEIN-PRIMAKOFF TRANSFORMATION 

For the sake of completeness we derive in this Ap- 
pendix explicit expressions for the different spin-spin cor- 
relators used in this work 



■/3 



(1) 



For this purpose, we make use of a Holstein-Primakoff 
transformation 



^25*^/1 - ^Oi, and 



(2) 



in the limit rii <^ 2S, with a,; satisfying bosonic com- 
mutation relations and rii = a|ai [ij. The creation op- 
erators al and annihilation operators ai satisfy bosonic 
commutation relations and the associated particles are 
called magnons. The corresponding Fourier transforms 



are straightforwardly defined as a, 



In harmonic approximation, the Heisenberg Hamiltonian 
Hp reads 



Hp 



(3) 



where Cq = uj^ + Ap ~ 4JS'[3 — {cos{qx) + cos{qy) + 
cos{qz))] + Ai? is the spectrum for a cubic lattice with 
lattice constant a = I and the gap Af is induced by the 
external magnetic field or anisotropy of the ferromagnet. 



TRANSVERSE CORRELATORS {S+{t)SZ^{0)) 

Let us now define the Fourier transforms in the har- 
monic approximation 

5Z„ = ^ y e'^^'^S- e^qr. ^ V25a_q . (4) 

From this it directly follows that 
C+-(t,q) = (5+(t)5Zq(0)) 



'q J 



(5) 



with eq w Dq^ + Aj? in the harmonic approximation. 



The Fourier transform is then simply given by 
C+- (c., q) = -= / dte-'-'C+- {t, q) 



/2tt 



■\/27r5(eq— a;) 



(6) 



The corresponding correlator in real space is then simply 
given by (q := |q|) 



^^"(-''■)=(^/dqe-^C+-(c.,q) 



(7) 



(^25^^ / dq^Pq^+AF-u.)e- 

- / dqdxqH{Dq^ + Ap - a;)e*«'' 
1 J~i Jo 



45 1 



/•OO 

— / dqq6{Dq^ + Ap - i^j) sin(gr) . 
1 Jo 



r e'^'^ - 1 Jo 
Let us now perform the following substitution 

V = Dq\ 

which gives for lo > Ap 
4S'/r 



(8) 



C+-(c.,r) 



2L»(eP 



1) 



dy5{y + Ap — 



(9) 



_ 25* 1 sin(^(w- AF)/£'r) 
^ "d e^" - 1 r 

We remark that 

C+"(cj,r) = 0, Lo<Ap. 



(10) 



We note the diverging behavior of the above correlation 
function for Ap = and w — > 0, namely 



lim 



1 sm ( 



e^^" - 1 



lim 



1 1 



(11) 



Similarly, it is now easy to calculate the corresponding 
commutators and anticommutators. Let us define 



5^(t,q) ■.^\{S+{t).SZ^m- 



(12) 



2 



It is then straightforward to show that 



(13) 



and therefore 



s 



\/27r J -co 



e»(^c,-c.)t(1^2nq) 
1 + 2. ' 



(14) 



FoUowing essentiaUy the same steps as the one performed 
above, we obtain the 3D real space anticommutator for 

I 

In real space, for the three-dimensional case, we obtain 



S'i°(w,q) = S'coth(/3w/2) X (15) 

S shU^/^^A^JjDr) 

= — cothfpa;/2) . 

D r 

(16) 

Let us now finally calculate the transverse susceptibility 
defined as 

X±(t,q)--*0(t)[5+(i),5lq(O)]. (17) 



As before, in the harmonic approximation, one finds 

In the frequency domain, we then have 
2iS 



(18) 



X-L('^,q) 



2tt Jo 
2S 



1 



(19) 



V27r eq - w + 177 
and thus in the small q expansion 

, , 25 1 
X-L(w,q) = 



2tt -Dq2 + Ap - to + iij 



(20) 



Xf(-,r) = - 



25 27r 



oo /•! 



2^ (27r)3/2 7o 
45 2tt 1 



dxdqq^ 



1 



2^ (27r)3/2 r 



1 

dqq 



Dq^ + Ai? — w + 177 
1 



Dq^ + Af - w + irj 



sa\{qr) . 



(21) 



Making use of the Plemelj formula we obtain for cj > Ap- 



Xf (c^,r) = - 



dqq 



Dq^ + Ap -UJ + bj 



sm{qr) 



2S 2-K I f°° , 1 
V2^ (2^)3/2 r 

25 27r 1 / - ^ q . , ^ • 

:-F I dq—^ sm(gr) + i , — ,^ ^ 

,Dq^ + AF-Lo ' ^/2^(27r) 



V2^ (2^)3/2 r 
S co& {r^{uj~ Af)/D) . 5 &in{^{uj ~ AF)/Dr) 



25 27r2 1 



2tt 1 f 
.3/2 - / dqqS{Dq'^ + Af - uj)sm{qr) 
) ^ J —00 



D 
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(22) 
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It is worth pointing out that the imaginary part of the where If — 



susceptibility vanishes 



(23) 



and therefore the susceptibility is purely real and takes 
the form of a Yukawa potential 



5 e^'^/^p 

Xf (c^,r) = --— — , u:<Af, 



(24) 



AiJ'— a; ' 

Note also that the imaginary part of the transverse sus- 
ceptibility satisfies the well-know fiuctuation dissipation 
theorem 



5f (c.,r)=coth(/3c./2)xf (c.,r)" 



(25) 



In three-dimension the susceptibility decay as l/r, 
where r is measured in lattice constants. For distances 



3 



of order of Ifim this leads to four orders of magnitude 
reduction. 

For quasi one-dimensional ferromagnets such a reduc- 
tion is absent and the transverse susceptibility reads 



S 



Xl°(c.,r) = -^;^e-'-/'-, w<Af, (26) 

where I p is defined as above and the imaginary part van- 
ishes as above, i.e., 



Similarly for lu > Ap we have 



sin \^^/{uj - Ap)/Dr 
^D{oj-Ap) 



(27) 



(28) 



and 



xl°(^,r)" = 77H 



S 



D 



2DV uj- A, 



■ cos 



V(w - Ap)/Dr^ 



(29) 



Applying Wick's theorem and performing a Fourier 
transform, we obtain the susceptibility in frequency do- 
main 



XII (w,q) = - 



N uj - eic+q + e-k + iri 



(31) 



where rik is the magnons occupation number, which is 
given by the Bose-Einstein distribution 



nk 



(32) 



LONGITUDINAL CORRELATORS (Sq(t)Siq(O)) 
The longitudinal susceptibility reads 

X|l(^,q) = -^^^(^)[5^(^),5lq(0)] (30) 



where ek is again the magnon spectrum (ek = Wk + A^? « 
_Dk^ + Ap for small k). Note that the longitudinal sus- 
ceptibility is proportional to 1/5, due to the fact that 

Ek — Ck+q = t^k — t^k+q OC S. 

Since we are interested in the decoherence processes 
caused by the longitudinal fluctuations, we calculate the 
imaginary part of X|| (i^, q) which is related to the fluctua- 
tions via the fluctuation- dissipation theorem. Performing 
a small q expansion and assuming without loss of gener- 
ality w > 0, we obtain 



q',q" 



xf^{^^ q)" = -^^^ j cfk(nk - nk+q)<5(a;k - Wk+q + t^) 
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(5 fc- 



2Bqx 



2Dqx 
1 



1 
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pi uj+Af+D{ 



(33) 



-1/ 



Next, since we are interested in the regime where uj ^ T more, we approximate the distribution function = 
(and thus /3tj > 1), we have > rife+g. Further- "'1'^^/^'^^'' (this is valid when Buj^ < 1) and arrive at 

the following expression 
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2Dqx 



bj - Dq^ 
2Dqx 



2 -^(a. + d(^)^) 



-Ei 



-0Ai 



2Dqx 



i ( -4 - PDq^ + 2l3uj - 



(34) 



where Ei(2;) is the exponential integral function. We also 
need the the real space represenation obtained after in- 
verse Fourier transformation, 
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(c.,r)" = / dggxr(^'9)"sin(gr). (35) 



In order to perform the above integral we note that the 
imaginary part of the longitudinal susceptibility, given by 
Eq. (|34p , is peaked around q = ^JZTpD with the width of 
the peak [l/y/pD] much smaller than its position in the 
regime we are working in {uj ^ T). For r = 0, the 
integration over q can be then performed approximately 
and yields the following expression 



2P'^D^ 

,/3Af _ 



(36) 



X Erfc(e^'^'^^/Ve'3^^ - 1)) 



where Erfc(z) denotes the complementary error function. 
It is readily observed from the above expression that the 
longitudinal fluctuations are exponentially suppressed by 
the gap. Assuming that Aj? ^ T, we obtain the following 
simplified expression 

^tF — e7rErfc(l) 



w, r 



0)'' 



(37) 
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We observe that, since J(w) = x\\{^^'^)" i tli6 longitudi- 
nal noise of ferromagnet is — as the transverse one — sub- 
ohmic 0- 

Next we calculate the longitudinal fluctuations for the 
case of a quasi- one-dimensional ferromagnet {/S.p 3> T) 
and obtain 



0)" = ^ 
47r 



dk 



1 



dq 
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6{uj - Dq^ - 2Dkq) 
1 



SCHRIEFFER- WOLFF TRANSFORMATION 
AND SPIN-SPIN COUPLING 

Exchange coupling to the ferromagnet 

The total Hamiltonian of the system reads 

H = Hf + H, + AY, (li'^^S- + a-S+) + , 

(39) 

where we identify the main part as = Hp + and 
the small perturbation as the exchange coupling V ~ 
A (Ti ■ Sn ■ The Hamiltonian of the ferromagnet reads 
Hp = —J^(^r Sr ■ Sr', while the Hamiltonian for the 

two distant qubits is H^ = ^ J2i=i 2 '^f ■ 

The second order effective Hamiltonian Q is given by 
H^l^ ^Ho + U, where 



U = — hm , 

2 v^o+ Jo 



dte-'^'[V{t),V] , 



(40) 



where Vit) = e'^^Ve 
We have 



1 - cos(Ai) 



(7, - i sin(At)CTf , 
(41) 



anda-(i) = a+(t)t. 

Recalling that the zz susceptibility can be neglected 
and that only the transverse susceptibility contributes, 
we obtain the following result from Eq. (pUj) . U = 

Y.^J^ dte-^'Y.^,U,j 



u,. 



[a-{t)Sl{t),a+S-\+\^.c 



(42) 



Finally, by rewriting cos(At) 



g/3(w+AF-i-r>fe2) _ 2 

I _ e-'^^^^H^flZJfc^ ]>yp^^^^y^ transverse spin susceptibility 

/3Af /oo\ 

X-l('^.i"j - Tj) = -i lim , 



-, sin(At) 



, and using the definition of the real space 



—00 
7 



where 7 is a numerical factor of order 1 . 



(43) 
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we obtain by inserting Eq. (|4T|) into Eq. (|42) 

f ^^jO) x±(A)+X±(-A) 
8 ^ V 2 4 



^ ij ^ ^ 

A2^x±(A)-x±(-A) 



■E 



erf 0-+ + h.c. 



(44) 



Since the decay length of the susceptibihty r) is large 
only close to the resonance, Ap A, we can simplify 
the above equation by neglecting x(~ A, r) and x(0, r) in 
comparison to x(A,r) which is assumed to be close to 
the resonance. Within this approximation we arrive at 
Eq. (3) of the main text. 



Dipolar coupling to the ferromagnet 

The total Hamiltonian of the system reads 

H = Hf + H^ + Hi, (45) 

where the interaction between the spin-qubit and the fer- 
romagnet, Hj, reads 



[^^■Sr Z2 ) ' (46) 



where /i is the magnetic moment of the ferromagnet, 
which is usually on the order of the Bohr magneton /i b ■ 
After performing a second order SW transformation we 
obtain the correction to the uncoupled Hamiltonian 



U : 



lini — 



1 



3(<Ti(i)T')(5r'(t)-r') 



(47) 



(T2 • Sr" — 



1 ^ 2, 



where a = [V/NY/'^ is the lattice constant of the ferro- 
magnet, This constant emerges in the above expression 
when going from the sum over the ferromagnetic spins to 
the integral over the ferromagnet volume. The above ex- 
pression can be significantly simplified if a specific geom- 
etry of the ferromagnet is assumed. Namely, we assume 
a dog-bone shaped ferromagnet — two ferromagnetic discs 
joined by a thin ferromagnetic channel. We assume that 
the radius of the two discs is much smaller than the dis- 
tance between their centers (i?o ^ L). Within this as- 
sumption we can take for the susceptibility between two 
points at opposite discs the same as the ID susceptibility 
given in Eq. (j29|) . Now there are two independent axial 
symmetries under the rotations around the axes passing 



through the spin qubit and the adjacent center of the 
disc (we assume the z-axis to be parallel to these axes). 
Due to this symmetry only terms like r'j^r'_ and r"r" 
will give non-zero results upon integration over the discs. 
We also neglect the longitudinal susceptibility compared 
to the transverse one and arrive finally at the following 
expression 



cff 
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xi°(A,L)(2afa| 



-'^1^2), (48) 



where we retained only the on-rcsonance susceptibility, 
i.e., we neglected Xj^(~A) and xi^(O) in comparison to 
Xj^(A). The constant B in the above formula reads 



B = 



2a3 



d2 



pdp 



1 



(z2+p2)3/2 



2 - 



3p2 



^2 _|_ p2 



(49) 

where d2 — di is the thickness of the discs, di is the dis- 
tance in z-direction from the spin qubit to the bottom 
of the adjacent disc. The above integral can be easily 
performed and yields the following result 



B 



2a3 



di 



V4 



Rl 



(50) 



IMPLEMENTATION OF TWO-QUBIT GATES 

Two qubits interacting via the ferromagnet evolve ac- 
cording to the Hamiltonian H = Hcs+Ha, the sum of the 
qubit-qubit coupling and Zceman term. These contribu- 
tions, by and large, do not commute, making it difficult 
to use the evolution to implement standard entangling 
gates. Nevertheless, since acts only in the subspace 
spanned by {|+, +) , |-, -)} and A > = A^x±{^)/'^, 
we can neglect the effect of i?cfr in this part of the space 
and approximate it by its projection in the space spanned 
by vectors {|+, -) , |-,+)} 



(51) 



Within this approximation, the coupling in i/^fj and Zee- 
man terms now commute. 

We consider the implementation of the iSWAP gate 
t^iSWAP = Q-ii'^i '^2 +'^1 '^2)3^/4^ which can be used to im- 
plement the CNOT gate. The Hamiltonian H' already 
contains a crfal^ term, so implementation of the iSWAP 
gate requires only that the erf erf term is transformed into 
erf erf by appropriate local rotations. This is achieved 
with the following sequence 

c/iswAP = nin2e'"''e~'"'^«'nin2, (52) 

where t = 37r/(4Ji2) and H denotes the single qubit 
Hadamard rotation. When iSWAP is available, the 
CNOT gate can be constructed in the standard way 



UcNOT 



-<^i e*4'^2 e^-<^2 C/iswApe-*-'"i C/iswApe'"'"^ . 

(53) 
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Since H'^g is an approximation of HaS, the above se- 
quence will yield approximate CNOT, J7cnotj when used 
with the full Hamiltonian. The success of the sequences 
therefore depends on the fidelity of the gates, ^(C^cnot)- 
Ideally this would be defined using a minimization over 
all possible states of two qubits. However, to character- 
ize the fidelity of an imperfect CNOT it is sufficient to 
consider the following four logical states of two qubits : 
\+, 0) , |+, 1) , | — , 0) , and | — , 1). These are product states 
which, when acted upon by a perfect CNOT, become the 
four maximally entangled Bell states j^"*") , \'^'^) , l^"") , 
and l^*^), respectively. As such, the fidelity of an imper- 
fect CNOT may be defined, 

FiUimT)^ , min | («, j| C/cnot^cnot K.?) 
«e{+,-} je{o,i} 

(54) 

The choice of basis used here ensures that ^'(?/cnot) 
gives a good characterization of the properties of C^qj^qt 
in comparison to a perfect CNOT, especially for the re- 
quired task of generating entanglement. For realistic pa- 
rameters, with the Zeeman terms two order of magni- 
tude stronger than the qubit-qubit coupling, the above 
sequence yields fidelity for the CNOT gate of 99.976%. 

To compare these values to the thresholds found in 
schemes for quantum computation, we must first note 
that imperfect CNOTs in these cases are usually mod- 
elled by the perfect implementation of the gate followed 
by depolarizing noise at a certain probability. It is known 
that such noisy CNOTs can be used for quantum compu- 
tation in the surface code if the depolarizing probability 
is less than 1.1% Q. This corresponds to a fidelity, ac- 
cording to the definition above, of 99.17%. The fidelities 
that may be achieved in the schemes proposed here arc 
well above this value and hence, though they do not cor- 
respond to the same noise model, we can expect these 
gates to be equally suitable for fault-tolerant quantum 
computation. 



DECOHERENCE 

The two-qubit system dynamics can be described by 
the reduced density matrix pK{t) = tr(p(t)), where p is 
the density matrix for the total system and the trace is 
taken over the degrees of freedom of the fcrromagnct. 
The reduced density matrix pBit) satisfies the General- 
ized Master Equation (GME) Q 

pnit) = [Ha,PR{t)] + f dt'J:{t - t')p^{t') . (55) 



The self-energy Yi{t — t') is a super-operator acting on 
PK{t). We assume that at the initial time t = Q the 
system of the qubit spins and the system of the ferro- 
magnet are decoupled, and the ferromagnet is in ther- 
mal equilibrium, p{Q) = /5r(0) ® p'^'^ , where p'^'^ ~ 



We calculate the self-energy Yi{t — t') of Eq. (1551) in 
Born approximation. We assume that the reduced den- 
sity matrix p^{t) varies slowly compared to the time scale 
of the ferromagnet, h/ J ^ and thus we employ the Markov 
approximation in Eq. (|55l) . The result of this reads 



Pk = -i[Ha,PR 




PR , 



(56) 



where Vp {p = 1,2) describes uncorrclatcd noise (i.e. Ko- 
rringa relaxation), and Pcorr describes correlated noise 
between the two qubits under the consideration. We ig- 
nore the correlated noise in what follows and set Pcorr = 
in Eq. ([55)1 . The reason for such a treatment is that we 
obtain the decoherence times due to uncorrelated noise to 
be much bigger than all the other time scales in the sys- 
tem, and the inclusion of the correlated noise cannot 
alter this statement Q. We have ignored the Stark and 
Lamb shift since these second order corrections are much 
smaller than the first order one we already included, see 
Eq. (2) of the main text. Therefore, in the main text we 
study decoherence due to coupling a single qubit to the 
ferromagnet. 

EXACT TREATMENT 

As noted earlier, the transversal noise of the ferromag- 
net is gapped and it does not contribute to decoherence 
in the lowest order. This motivate us to study a model 
where we have only longitudinal coupling to the ferro- 
magnet. After making use of harmonic approximation 
for the ferromagnet we arrive at the Hamiltonian that 
can be exactly solved since it can be related to the ex- 
actly solvable 'independent Boson problem' Instead 
of Eq. (1) of the main text, we have a simplified Hamil- 
tonian of the system 



H = Ho 



v. 



(57) 



with i^O = Eq ^qaqftq + £<^z and V = Eq ^(Oq + flq) ■ 

Since CTj commutes with the total Hamiltonian, only 
dephasing can occur in the system. In order to study 
dcphasing we have to calculate the following quantity [10[ 



(a_(t))=e-*/^a_(0))x 



(58) 



Tcxp ( - 



Vdt' I Texp ( - 



Vdt' 



-Hp/knT 



/tr (e" 



^HF/kBT\ 



The average in the above formula can be calculated using 
a cluster expansion |ll| . and since the perturbation V is 
linear in the bosonic operators, only the second order 
cluster contributes. Therefore, the final exact result for 
the time-evolution of cr_ (t) reads 
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where S{t) = {[V{t),V{0)]+). Furthermore, we can ex- 
press this result in terms of the spectral function J{io) — 
tanh(;3w/2)5(w) 



dephasing is much stronger then in the case when the 
splitting is present. This was the reason for including 
the splitting of the qubit. 



-it)) 



X exp 



(0))x 

1 f duj 
2^ 



J(a;)coth(^cj/2) 



sm 



(60) 
M/2) \ 

(c./2)2 ; • 



Note that the above formula is of the exactly same form 



as the one for a classical Gaussian noise [12| . Now we 
apply this formula to the case of the ferromagnetic bath. 

At zero temperature, we obtain from Eq. ((60)) the de- 
phasing behavior described by the evolution {a-{t)) ^ 

g " , In stark contrast to this, a naive use of the 

Born-Markov approximation leads to the incorrect con- 
clusion that the dephasing is absent Q. However, we note 
that the latter approximation is no longer self-consistent 
in this regime since the perturbation (coupling to the fer- 
romagnetic bath) becomes bigger than the system energy 
scales (which vanishes for zero splitting). Since the lon- 
gitudinal noise of the ferromagnet is completely absent 
at T = 0, we next analyse dephasing for non-zero tem- 
perature. The Born-Markov approximation again leads 
to divergent result at non-zero temperature, which again 
is spurious. The exact solution for T 7^ 0, and for long 
times t ^ h/T leads to the well-behaved dynamics of the 
form 



(61) 
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